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We investigate the universality of multi-spin systems in architectures of various symmetries of 
coupling type and topology. Explicit reachability sets under symmetry constraints are provided. 
Thus for a given (possibly symmetric) experimental coupling architecture several decision problems 
can be solved in a unified way: (i) can a target Hamiltonian be simulated? (ii) can a target gate 
be synthesised? (iii) to which extent is the system observable by a given set of detection operators? 
and, as a special case of the latter, (iv) can an underlying system Hamiltonian be identified with a 
given set of detection operators? Finally, in turn, lack of symmetry provides a convenient necessary 
condition for full controllability. Though often easier to assess than the well-established Lie-algebra 
rank condition, this is not sufficient unless the candidate dynamic simple Lie algebra can be pre- 
identified uniquely, which is fortunately less complicated than expected. 
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Introduction 

Controlling the quantum dynamics of experimentally 
manageable systems is paramount to exploiting the great 
potential quantum systems inherently promise: they may 
provide access to performing computational tasks or to 
simulating the behaviour of other quantum systems that 
are beyond experimental handling themselves. Moreover, 
quantum systems can efficiently simulate classical sys- 
tems [l|, Q sometimes separating controllable parame- 
ters in the quantum analogue that classically cannot be 
tuned independently. — Both in simulation and compu- 
tation the particular advantage becomes obvious when 
the complexity of a problem reduces upon going from a 
classical to a quantum setting Q . On the computational 
end, most prominently, there is the exponential speed- 
up by Shor's quantum algorithm of prime factorisation 
[1, 0| relating to the class of quantum algorithms 0, 0| 
efficiently solving hidden subgroup problems While 
the demands for accuracy ('error-correction threshold') 
in quantum computation may seem daunting at the mo- 
ment, the quantum simulation end is by far less sensi- 
tive. Thus simulatin g qu antum systems [3| — in particular 
at p hase-transitions — has recently shifted into focus 
[ill . [l^ . [isl . [l3 | mainly inspired by promising experimen- 
tal prog ress in cold atoms in optical lattice potentials 
[H, ^ as well as in trapped ions In order 

to assess feasibility, Kraus et al. have explored whether 
target quantum systems can be universally simulated on 
translationally invariant lattices of bosonic, fermionic, 
and spin systems [l^. Their work can also be seen as 
a follow-up on a study by Schirmer et al. [131 specifically 
addressing the question of controllability in systems with 
degenerate transition frequencies. 

Quite generally, quantum control has been recognised 
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as a key generic tool [IH needed for advances in exper- 
imentally exploiting quantum systems for simulation or 
computation (and even more so in future quantum tech- 
nology). It paves the way for constructively optimising 
strategies for experimental implemention in realistic set- 
tings. Moreover, since such realistic quantum systems are 
mostly beyond analytical tractability, numerical meth- 
ods are often indispensable. To this end, gradient flows 
can be implemented on the control amplitudes thus it- 
erating an initial guess into an optimised pulse scheme 
[H; l23[ [23l . This approach has proven useful in spin sys- 
tems [25fas well as in solid-state systems Moreover, 
it has recently been generalised from closed systems to 
open ones [23|, where the Markovian setting can also be 
used as embedding of explicitly non-Markovian subsys- 
tems [111. 

Clearly, for universal quantum computation or simu- 
lation, the quantum hardware (characterised by the sys- 
tem Hamiltonian) and the quantum software (operating 
via the interface of control Hamiltonians) have to com- 
bine so that any unitary target operation can be per- 
formed irrespective of the initial state of the quantum 
system. More precisely, this means the system has to 
be fully controllable (as has been pointed out in differ- 
ent contexts [29l. [30|). Yet often the quantum hardware 
comes in designs with a certain symmetry pattern re- 
flecting the experimental set up. For instance, this is 
the case in quantum lattices, in quantum networks, or in 
spin chains serving as quantum wires [3l| for distributed 
quantum computation js^]. However, symmetric pat- 
terns may have crucial shortcomings since symmetry may 
preclude full controllability. On the other hand, avoid- 
ing symmetry-restricted controllability need not be com- 
plicated from an experimental point of view: in Ising 
spin chains it soon emerged that polymers {ABC)n made 
of three different qubit units A,B,C are fully control- 
lable [H, [Hi. Later irregular ABAAA... systems of 
just two qubit types A, B [33| and even A — A - ■ ■ A — B 
systems turned out to be fully controllable as well. In 



2 



these systems, spin qubits are meant to be locally con- 
trolled by operations that act jointly on all qubits labelled 
with the same letter and independently from controls on 
qubits with a different letter. More specifically, quantum 
systems coupled by Heisenberg XXX type interactions 
turned out to be fully controllable when the local con- 
trols are Hmited to a small subset of qubits, as with time 
these actions can then be 'swapped' to neighbouring spins 

[m, M, M, m, m, so, SI. 

Given the power of Lie theory to assess controllability 
as well as observability in classical systems [H, S3|, such 
a gradual case-by-case development asks for a more sys- 
tematic investigation on the quantum side. Here we ad- 
dress controllability and observability with and without 
symmetry restrictions in multi-qubit systems, where the 
coupling topology is generalised to any connected graph 
[44 . .45, . .461] . In doing so we also specify reachability sets 
as subgroup orbits for systems that are not fully control- 
lable since it is important to know for which dedicated 
tasks of simulation or computation symmetry restricted 
systems can still be used. 

Scope and Overview: On a more general scale, it is 
important to be able to separate questions of existence 
(e.g.: is the system fully controllable? is it observable? 
can it be used for universal Hamiltonian simulation?) 
from questions of actual implementation {how does one 
have to steer a given experimental setup for implementing 
a target task with highest precision ?) . Otherwise the ele- 
gance of constructive proofs of existence may all too often 
come at the cost of highly 'suboptimal' implementation 
when translated directly into experimental schemes ir- 
respective of the actual given setting. Therefore, here 
we advocate to exploit the power of Lie theory for a 
unified framework addressing controllability (and within 
the context of tomography also observability) in a first 
step, while resorting to quantum control in a second step 
for actual implementation optimised for the given exper- 
imental set up. 

In view of practical applications, we also systematically 
explore what we call task controllability, i.e. which tasks 
are feasible on which type of quantum system hardware. 
In turn, the unified picture also provides the rules how 
to design quantum hardware for a specific task. 



I. CONTROLLABILITY 

Here we address Markovian dynamics of quantum sys- 
tems, the free evolution of which is governed by a system 
Hamiltonian Hd and, in the case of open systems, by 
an additional Markovian relaxation term T taking GKS- 
Lindblad form. Whenever we talk about controllability, 
we mean full operator controllability thus neglecting more 
specialised notions like pure-state controllability [47[. 



A. Quantum Dynamical Control Systems 

The interplay between the quantum system and the 
experimenter is included by control Hamiltonians Hj ex- 
pressing external manipulations in terms of the quan- 
tum system itself, where each control Hamiltonian can be 
steered in time by control amplitudes Uj{t). With these 
stipulations the usual equations of motion for controlled 
quantum dynamics can be brought into a common form 
as will be shown next. 

To this end, consider the Schrodinger equations 

m 

m)) = -t{Hd + Y,Mt)Hj)\m) (1) 

m 

U{t) = -i{Hd + Y,u,{t)H,) U{t) , (2) 

where the second identity can be envisaged as the oper- 
ator equation to the first one. It governs the evolution 
of a unitary map of an entire basis set of vectors rep- 
resenting pure states. Using the short-hand notations 
Hu := Ha + J2]Li^A^)^j and ad//(-) [H,{-)] also 
consider the master equations 

m = -i[H„,p(t)]-r(p(t)) = -(iadH„ +r)p(t)(3) 

F{t) = -(^adH„ +r) o F{t) . (4) 

While p € i)C'c{N), F denotes a quantum map in GL{N'^) 
as linear image over all basis states of the Liouville space 
representing the open system. 

All these equations of motion ressemble the form of a 
standard bilinear control system (S) known in classical 
systems and control theory reading 

m 

X{t) = {A + J2nj{t)B,) X{t) (5) 
i=i 

with 'state' X{t) e C^, drift A e MatAr(C), controls 
Bj £ Mat AT (C ) , and control amplitudes Uj & R . For 
simplicity consider for the moment its linear counterpart 
with v&C^ 

X{t) ^ AX{t) + Bv (6) 

which is known to be fully controllable if it obeys the 
celebrated rank condition (see, e.g., (ist p 

raTik[B, AB, A^B,. A^-^B]= N . (7) 

Now liftin g th e bihnear control system (S) to group man- 
ifolds 0,1531 by X{t) e GL{N,C ) under the action of 
some compact connected Lie group K with Lie algebra ? 
(while keeping A, Bj G Mat Ar(C )), the condition for full 
controllability turns into its analogue known as the Lie 
algebra rank condition (i^ . S3, 

(A,Bj|.7 = l,2,...,m)Lic = « , (8) 

where (•)Lic denotes the Lie closure obtained by repeat- 
edly taking mutual commutator brackets. 
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Algorithm 1 [52] : Determine Lie closure for n-qubit system 
with given set of drift (or system) and control Hamiltonians 

1. Start with the inital basis B,^ := {Hd', Hi, ■ ■ ■ , Hm}- 

2. If m + 1 = dim su(2") ^ terminate. 

3. Perform all Lie brackets Ki = [Hj,Hk]. 

4. For each new Ki check linear independence 
from spanS^. If nothing new found terminate. 

5. Extend basis by all independent new {Ki} 
B.+i - {{Ki},Hi}. Goto (1). 




Figure 1: General coupling topology represented by a con- 
nected graph. The vertices denote the spin-i qubits, while 
the edges represent pairwise couplings of Heisenberg or Ising 
type. Qubits of the same colour and letter are taken to be 
controlled by joint local unitary operations as in Tab. I (or 
none: see Tab. II), while qubits of different kind can be con- 
trolled independently. For a system to show a certain sym- 
metry (brought about by permuting subsets of qubits), it is 
necessary that both the graph as well as the system plus all 
control Hamiltonians are invariant, see text. 

columnwise as vec (Ki) [53|. Our results have been cross- 
checked with GAP 4.4.12 |5J|]- 

Notation: Coupling Graphs 

Here we adopt the notation to represent a system by 
a graph G{V, E) whose vertices represent the qubits and 
whose edges denote pairwise coupHngs of Ising or Heisen- 
berg nature. In the subsequent examples, the qubits 
are taken to be locally controlled by operations that act 
jointly on all qubits labelled with the same letter and in- 
dependently from controls on qubits with a different let- 
ter. We assume that every qubit type is fully controllable 
so all SU (2)-actions are jointly admissible on each qubit 
with the same letter. Let us emphasise that although the 
current set of examples is taken with some local control 
on each qubit, the algebraic approach put forward here 
extends to systems with subsets of uncontrolled qubits 
coupled by isotropic Heisenberg interactions thus repro- 
ducing all the results on transmitting controls through 
infective graphs given recently [43 |. 

Characterisation by Symmetry 

In the following, we will characterise systems of re- 
stricted controllability in terms of symmetries. In the 
present setting, a Hamiltonian quantum system is said to 
have a symmetry expressed by the skew-Hermitian oper- 
ator s e 5u(A^), if 

[s,iJ^]=0 for all t/G {d;l,2,...,m} . (11) 

More precisely, we use the term outer symmetry if s gen- 
erates a SWAP operation permuting a subset of spin 
qubits such that the coupling graph is left invariant. In 
contrast, an inner symmetry relates to elements s not 
generating a SWAP operation in the symmetric group of 
all permutations of spin qubits in the system. 



B. Closed Quantum Systems: Full Controllability 
and Symmetry- Restricted Controllability 

Clearly in the dynamics of closed quantum systems, 
the system Hamiltonian Hd is the only drift term, 
whereas the Hj are again the control Hamiltonians. To 
fix notations, in systems of n qubits define N := 2", 
so these traceless spin Hamiltonians i Hy £ su{N) each 
generate a one-parameter unitary group of time evolution 
U^it) := e-''"" e SU{N). 

Transferring the classical result [i^l to the quantum 
domain j5l|, the bilinear system of Eqn. [2] is operator 
controllable if and only if the drift and controls are a 
generating set of su(A^) 

(iifd,«Hj|.7 = l,2,...,?7i)Lic = « = su(7V) . (9) 

In fully controllable systems to every initial state po 
the reachable set is the entire unitary orbit C'u(po) '■= 
{UpoW I U e SU{N)}. With density operators being 
Hermitian this means any final state p{t) can be reached 
from any initial state po as long as both of them share 
the same spectrum of eigenvalues. 

In contrast, in systems with restricted controllability 
the Hamiltonians only generate a proper subalgebra 

{iHd,iHj\j ^l,2,...,m)uo^tC5n{N) . (10) 



Algorithm for Computing the Lie Closure 

Given an n-qubit bilinear control system characterised 
by the drift and control Hamiltonians {Hd] Hi, ... , Hm}- 
Then the canonical Algorithm 1 tabulated above con- 
structively determines a basis of the associated dynamic 
Lie algebra ? [s^l- — Our implementation codes the or- 
thogonal basis of tensor products of Pauli matrices as 
quaternions simply represented by the Clifford algebra 
C^2(R) of quarternary numbers {0, 1,2,3} plus the Clif- 
ford multiplication rule allowing to calculate Lie brackets 
without matrix representations. The only time consum- 
ing step is the final rank determination by QR decompo- 
sition of a sparse coefficient matrix C £ Mat 4-7 collecting 
all the expansion coefficients to the Ki of Algorithm 1 
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In either case, a symmetry operator is an element of 
the centrahser (or synonymously the commutant) 

{H,y := {s G 5u{N) I [s, H,] = Vi^ € {d; 1, 2, . . . , m}} . 

Recall that the centraliser or commutant of a fixed subset 
m C g with respect to a Lie algebra g consists of all 
elements in g that commute with all elements in m. By 
Jacobi's identity [[a, 6],c] + [[6, c],a] + [[c, a],6] = one 
gets two properties of the commutant pertinent for our 
context: First, an element s that commutes with the 
Hamiltonians {iH,y} also commutes with their Lie closure 
i. For the dynamic Lie algebra t we have 

r := {s e 5u{N) I [s, k]^0 Vfc e i} (12) 

and {iHt,}' ~ t' . Thus in practice it is (most) convenient 
to just evaluate the commutant for a (minimal) generat- 
ing set {iH^} of t. Second, for a fixed k G i an analogous 
argument gives 

[si,fc]=0 and [s2,k]=0 [[si,S2],k] , 

so the commutant i' forms itself a Lie subalgebra to 
su{N) consisting of all symmetry operators. Since the 
commutant i' to the dynamic Lie algebra t is invariant 
under t, it is a normal subalgebra or an ideal to ? 

% i'] c t' . (13) 

NB: while the embedding Hilbert space Ti. of Hermitian 
operators is endowed with a weak (and strong) operator 
topology thus giving von Neumann's bicommutant the- 
orem [5a| H" = H well-known in C* -algebras [H^ . [57| . 
no such theorem extends to Lie algebras: there may be 
several inequivalent Lie algebras sharing the same com- 
mutant, in particular when the commutant is trivial. 

Within the commutant one may choose its centre 
3 = t' n t" of mutually commuting symmetry operators 
then allowing for a block-diagonal representation in the 
eigenspaces associated to the eigenvalues (Ai, A2, . . . ,X() 
to {zi, Z2, . . . , z^} = 3. A convenient set of symmetry op- 
erators embracing the outer symmetries are the ones gen- 
erating SWAP transpositions of qubits: they correspond 
to the 52 symmetry and come with the eigenvalues -f 1 
(gerade) and —1 {ungerade) . A block-diagonal represen- 
tation results if all the SWAP transpositions that can be 
performed independently are taken as one entry each in 
the £-tuple (Ai, A2, . . . , Af), while all those that have to 
be performed jointly are multiplied together to make one 
single entry in the tuple. This procedure is illustrated 
below in Examples 1 and 2. 

Notation: Clebsch-Gordan Decomposition 

Observe that in our notation a block- diagonal repre- 
sentation of a Lie algebra 6 = su(A^i) © 5u(A^2) gener- 
ates a group K = SU{Ni) © SU{N2) in the sense of a 
block-diagonal Clebsch-Gordan decomposition, while the 




Figure 2: Ising spin chain with joint S2 symmetry. 
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Figure 3: The drift and control Hamiltonians of Ex. 1 take 
block diagonal form corresponding to the Ag and Au repre- 
sentation of the 52-symmetry group. 

abstract direct sum of the algebras has a matrix repre- 
sentation as the Kronecker sum g = su(A^i) © su(A^2) ■— 
su(A^i) ® Iatj -|- (8)Su(A^2) and generates a group iso- 
morphic to the tensor product G = SU(Ni) ® SU{N2). 
Note that the direct sum of two algebras g and () (each 
given in an irreducible representation) has itself an irre- 
ducible representation as Kronecker sum g © f) [s^ . 



Examples with Symmetry-Restricted Controllability 

Example 1: Joint S2 Symmetry 

First, consider Ising n-spin-^ chains with odd num- 
bers of spins such as, e.g., the one in Fig. [21 which has 
an Ising coupHng graph L5. Denote the centre spin by 
C. An S'2-symmetry leaves the coupling graph and thus 
all Hamiltonians of drifts and controls Hi, ... , Hm} 
invariant under the joint permutation j ^ {N — j -I- 1) 
V j = 1 . . . [r7,/2j . In the 5'2-symmetry-adapted basis, the 
Hamiltonians take block-diagonal form with two blocks 
of different parity, gerade and ungerade. The dimensions 
of the respective blocks are 

(m+a)/2 , ^ 

d, := 2" L,™ J3'— - (14) 

k=a ^ ^ 

du := 2"~dg (15) 

where m denotes the number of symmetric spin pairs and 
a := TO(mod 2). Moreover, each block represents a fully 
controllable logical subsystem. 

In the system depicted in Fig. [21 spins A and B can be 
jointly exchanged with spins A' and B' , since the system 
has a mirror symmetry with C in the mirror axis. As 
shown in Fig. [31 the S'2-symmetry adapted basis entails 
block diagonal representations of the Hamiltonians with 
one block being of size 20 x 20 and one block of size 
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Figure 4: Coupling topology allowing for individually inde- 
pendent permutation symmetries 11^^' and Hbb', which to- 
gether reproduce the 5'2-symmetry of Ex. 1. 
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Figure 5: The drift and control Hamiltonians of Ex. 2 take 
block diagonal form corresponding to the gg, uu, gu and ug 
parities of the individual permutations Haa' and Hbb'- 
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12 X 12. Note that all the symmetrised Hamiltonians are 
traceless within each of the two blocks. The Lie closure 
gives a Clebsch-Gordan (CG) decomposed dynamic Lie 
algebra 

Si =su(20)®su(12) so Ki S'[/(20)©5[/(12) (16) 
and thus the Lie rank is 542 = 399 + 143. 

Example 2: Individual Permutation Symmetry 

By introducing further Ising couplings between spins 
A and B' and spins A' and B, the L5 system of Ex. 1 can 
be turned into the one represented by a non-planar graph 
in Fig. m note that now the spin pairs A, A' and B, B' 
can be permuted individually. Thus the block-diagonal 
representation consists of four blocks corresponding to 
the parities gg, uu and gu, ug with the respective sizes 
18 X 18, 2 X 2 and twice 6x6, see Fig. O The indis- 
tinguishable pair AA' (and BB') can be looked upon as 
a pseudo spin-1 and pseudo spin-0 system, because sym- 
metrisation allows for adding their spin angular momenta 
in the usual Clebsch-Gordan way. The Lie rank of each 
pair is 4 instead of 15 in a fully controllable spin pair. 

Now the Lie rank of the total system is 364 = 323 -I- 
35 -I- 6, i.e. the sum of all blockwise Lie ranks — with two 
exceptions: (i) the second 6x6 block is not of full rank 
since it does not contain any coupling terms; instead, it 
only collects elements from two independent su(2) subal- 
gebras which arise from the AA' pair and the C spin thus 
giving a Lie rank of 6; (ii) the 2x2 block does not con- 
tribute since it reduplicates one of the independent su(2) 
algebras already occuring in the second 6x6 block, which 
becomes obvious as the matrix elements in both blocks 
only occur jointly. Finally, since in all Hamiltonian com- 
ponents all the blocks are independently traceless, we 
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Figure 6: Coupled qubit systems with subsets of uncontrolled 
spins that are coupled via a connected coupling topology of 
Heisenberg interactions. Their dynamic Lie algebras are listed 
in Tab. II. Note that the graph of (e) is 'non infective' [ifil ]. 
yet the system is fully controllable. 



have a CG-decomposed dynamic algebra 

li = su(18) © su (2) © su(6) © (su(2)j.^i © su(2)) (17) 
I I 

generating a dynamic group 

K2 = SU{\%)®SU (2) © S'[/(6) © (S'[/(2)j-^i ® SU{2)) , 
I I 

(18) 

where the index j=l denotes the spin-1 representation of 
su(2) C su(3) and the bracket connects the su(2) of spin 
C occuring in two copies {vide supra). 
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Table I: Dynamic System Algebras to Qubit Systems with Joint Local Controls on Each Type (Letter) and Various Couplings 



System 


Coupling Types Lie Rank Block Sizes Lie Ranks 


System Lie Algebra 


Trace = 




blockwise 




blockwise 



A 


-A 








ZZ, XX, XY 


9 


3, 1 


9, 1 


5(u(3) ^ 


33u(l)) 


no 


A 


-A 








XXX 


4 


3, 1 


4,1 


s(u(2) < 


33u(l)) 


no 


A 


-B 








-all- 


15 


4 


15 


su 


(4) 


yes 


A 


1 B 








-none- 


6 


4 


6 


su(2) { 


Bsu(2) 


yes 


A 


- B 


-A 






-all- 


38 


6, 2 


35, 3 


su(6) e 


Bsu(2) 


yes 


A 


- B 




~ A 




-all- 


135 


10, 6 


100, 36 


s(u(10) 


ffiu(6)) 


no 


A 

1 


- B 


- B 


-A 




ZZ, XX, XY 


135 


10, 6 


100, 36 


s(u(10) 


ffiu(6)) 


no 


A 
1 


— is 


— is 


A 

— A 




V v v 


115 


10, 6 


100, 16 


s(u(10) 


®u(4)) 


no 


A 


-B 


-c 


- B 


-A 


-all- 


542 


20, 12 


399, 143 


su(20) C 


Bsu(12) 


yes 


A 

1 


-B 


-c 


- B 


-A 


-all- 


543 


20, 12 


400, 144 


s(u(20) < 


©u(12)) 


no 


u 


-B 

1 


-c 


1 

- B 


-A 


-all- 


364 


18, 2, 6, 6 


323, (3), 35, 6 


given in lengthy Ean.l|17p 


yes 


A 

1 


-B 


-c 


- B 


-A 


-all- 


1023 


32 


1023 


su(32) 


yes 


A 


-B 


-c 


- C 


- B - 


- A all 


2079 


36, 28 


1296, 784 


s(u(36) 1 


©u(28)) 


no 



Table II: Dynamic System Algebras to Qubit Systems with Proper Subsets of Uncontrolled Qubits 



System 
shown in 


Coupling Types 
non-zero 3i 


Lie Rank 


Block Sizes 


Lie Ranks 
blockwise 


System Lie Algebra 


Trace = 
blockwise 


Pig. Wfi) 


^arb. 


15 


4 


15 


su(4) 


yes 


Fig. Wh) 




30 


4, 4 


15, 15 


su(4) esu(4) 


yes 


Fig. ms) 


XXX 1 ^ zz 


30 


4, 4 


15, 15 


su(4) su(4) 


yes 


Fig. lad) 


XXX 


255 


16 


255 


BU(16) 


yes 


Fig. He) 


Tl j2 , t3 

'JxXXt '-' XXX / ^ XXX 


255 


16 


255 


su(16) 


yes 


Fig.Hf) 


" 3^a:a: 5 zz 


126 


8, 8 


63, 63 


su(8) esu(8) 


yes 


Fig.Hg) 


7-1,2,3 t4 

XXX ! ZZ 


510 


16, 16 


255, 255 


su(16) esu(16) 


yes 



C. Task Controllability 



Clearly a set of non-trivial symmetry operators pre- 
cludes full controllability since the symmetry restrictions 
entail that the dynamic Lie algebra t is but a proper 
subalgebra of su(iV) . Thus 



(cxpt) K C S'C/lAf) (19) 



and the reachable sets take the form of subgroup orbits 

Ok{pq) := {Kpd<-^\ KeK] . (20) 

This is one reason why in the case of symmetry-reduced 
controllability it is advantageous to have specified the 
dynamic Lie algebra t explicitly by the methods intro- 
duced above. Another reason is to be able to assess 
the feasibility of specific tasks in the case of non uni- 
versal systems of reduced controllability. More precisely, 
a Hamiltonian quantum system characterised by {iH^} 
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is called to be task controllable with respect to a tar- 
get unitary module (or gate) Ut if there is at least one 
Hamiltonian Ht on some branch of the 'logarithm' of 
Ut so that Ut = e**^ • e~*^^ (with arbitrary phase (p) 
and iHT G t = {iH^)uc- For a detailed discussion of 
the phase and its experimental significance in systems 
with time hierarchies between local and coupling opera- 
tions see Ref. Needless to say only fully controllable 
systems are task controllable with respect to all unitary 
target modules. 



D. Lack of Symmetry versus Full Controllability 

Ultimately the question arises under which conditions 
the absence of any symmetry in turn implies full con- 
trollability. In the special case of pure-state controllabil- 
ity, this interrelation was analysed in In the gen- 
eralised context of full operator controllability the issue 
has been raised in [5^, inter alia following the lines of 
[60| . however, without a full answer. In particular, here 
we ask: what is special about quantum systems where 
the drift Hamiltonian comprises Ising or Heisenberg-type 
couplings in a topology that can take the form of any 
connected graph? — To this end, observe the following: 

Lemma I.l Let t C su(iV) be a matrix Lie subalgebra to 
the compact simple Lie algebra of special unitaries su{N) . 
If its commutant ( centraliser) t' of I in su{N) is trivial, 
then 

(1) t is irreducible; 

(2) t is simple or semi-simple. 

Proof. (1) The unitary representation as matrix Lie al- 
gebra t C su{N) is fully reducible by the Schur-Weyl 
theorem. Now as there is no invariant subspace other 
than the trivial ones, the representation is irreducible. 
(2) Since t is by construction a Lie subalgebra to the 
compact Lie algebra su{N) it is compact itself. By com- 
pactness it has a decomposition into its centre and a semi- 
simple part t = 3t ©s (see, e.g., [6l| Corollary IV. 4. 25). 
As the centre = n t is trivial, t can only be semi- 
simple or simple. ■ 

Moreover, by considering invariant subalgebras (ideals) 
one gets further structural properties. 

Lemma 1.2 Let the Lie closure t of a set of drift and 
control Hamiltonians {iH^} be a compact Lie algebra 
with trivial commutant t' in 5u{N). Then t is a simple 
Lie algebra, if there is a drift Hamiltonian iHd G {iHi,} 
satisfying the two conditions 

1. the coupling topology to iHd takes the form of a 
connected graph extending over all the qubits; 

2. the coupling type is Ising ZZ or Heisenberg XX , 
XY, XXX, XXY, XYZ. 



Proof. In view of applications in physics, we address two 
cases: (1) Consider the scenario where we have full lo- 
cal controllability on each type of qubit and a connected 
coupling topology of pairwise Ising (or Heisenberg) inter- 
actions (cp. Tab. I). With i' being trivial, the Lie closure 
of the controls without the drift {iH^} \ iHd is a semi- 
simple Lie algebra given in the irreducible representation 
t) = [)i ffi f}2 © ■ ■ ■ © [)r Si where the sum typically ex- 
tends over the local qubit spaces. As the drift term shows 
a coupling topology represented by a connected graph, 
there is no subalgebra f)fc that remains invariant (in the 
sense of being normalised) under all the coupling terms in 
iHd- Thus the decomposition into a Kronecker sum that 
exists for any semi-simple Lie algebra disintegrates when 
taking into account all the coupling terms and hence the 
Lie closure must turn into a simple Lie algebra at the 
latest upon including iHd- 

(2) Consider the scenario where we have full local con- 
trol on a (non empty) subset of qubits while others have 
no local control at all (cp. Tab. II). If one has again a con- 
nected coupling topology such that all the uncontrolled 
qubits have at least one coupling of Heisenberg XXX (or 
XXY or XYZ) type, while the coupling between pairs of 
locally controlled qubits may take Ising form as in (1), 
then observe the following: controlled local qubits trans- 
fer to their uncontrolled neighbours by the Heisenberg 
XXX (XXY, XYZ) interaction. Note that a commuta- 
tor of a single-qubit operator with H^xx acts like a SWAP, 
but the Heisenberg XXX interaction is entangling (cp the 
VSWAP) and thus does not allow for a Kronecker-sum 
structure giving a merely semi-simple system algebra; it 
enforces simplicity. Moreover if there is no symmetry 
in the coupling Hamiltonians, even several uncontrolled 
qubits coupled to one single controlled qubit maintain ir- 
reducibility (by a trivial commutant resulting from lack 
of symmetry) and simultaneously simplicity by coupling 
local qubit spaces (see Fig.[6](e) and Tab. II). This gen- 
eralises the scenario of 'infecting' graphs fi^l. ■ 

Example: Consider a two-qubit system with local con- 
trollability, so the closure to {iH^}\iHd is su(2) ©su(2), 
which is semi-simple (and isomorphic to so(4)) and has 
just a trivial commutant. Upon including, e.g., an Ising 
coupling, the Lie closure of the full {iH^} turns into 
su(4), which is simple. 

So lack of symmetry gives a trivial commutant, which 
in turn implies irreducibility and together with compact- 
ness it entails (at least) semi-simplicity, while a connected 
coupling topology on top finally ensures simplicity. 

Corollary I.l Let a controlled closed n spin-^ qubit 
quantum system be given which is characterised by the 
respective drift and control Hamiltonians {iHi,} C 5u(A'') 
and their Lie closure {iHd, iHj \i — 1, 2, . . . , m)Lie = ? 
being a subalgebra to 5u{N), where N > 4. 

If a check on the {H^} ensures the only candidates for 
a simple real compact Lie algebra i C su{N) are of the 
type ai (su(€ -I- 1)) in n spin-^ representation, then the 
following assertions are equivalent: 
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(a) the system is fully controllable; 

(h) the Lie closure is a representation t su(iV); 

(c) the Lie closure t is a simple Lie algebra and its 
commutant is trivial 

= {A • In} with Xei-R ; (21) 

(d) t is an irreducible representation of a simple Lie 
algebra in i ■ f)ero(A^). 

Proof. For (a) (6) see while (c) (d) is 
Schur's Lemma applied to the representation of Hilbert 
spaces, see, e.g.. Proposition 2.3.8. in [H^l- The forward 
implication (6) (c) <;=^ (d) readily follows since 

su{N) is simple and the Lie closure t = su(A^) leaves 
no non-trivial invariant subspace or subalgebra of skew- 
Hermitian operators in i ■ f)er(iV), hence it is irreducible. 

Only the backward implication (c) <;=4> (d) =4> (6) is 
intricate because su(A^) is not unique among the simple 
compact Lie algebras: they comprise the classical types 
aj.bi, C(, Of as well as the exceptional ones e(,fi and ge 
[63 |. all of which have to be ruled out along the lines of 
the subsequent constructive candidate filter just leaving 
the type ai. 

Lemma 1.3 (Candidate Filter) Given a set of Hamil- 
tonians {iH^} generating the Lie closure I C 5u(2") with 
the promise that t is an irreducible n spin-^ represen- 
tation of a compact simple Lie subalgebra to su(2") and 
n > 2. Then t = su(2") unless one of the following three 
conditions is satisfied: 

1. all the {iH^} are jointly conjugate to a set {iHy\, 
where each element is real and skew- symmetric 
Hi - -i?,, or 

2. all the {iH^} are jointly conjugate to a set {iH^}, 
where each element is unitary and symplectic 
JHy = —HlJ with J :— iuy ® 1^/2 so ~ — Iw, 
or 

3. the {iH^} generate an irreducible unitary subalge- 
bra 5u{N') C 5u{N) with N' < N that is compatible 
with an irreducible n spin-^ representation. 

Proof of the Lemma. For multi-spin systems with 
> 4 consider the compact classical simple Lie algebras 
[g^ I as candidate subalgebras to su(A^) 

ai > 1) : su(£ + 1,C) 

ii>2) : so(2£+ 1,R) 

ci{e>3) : sp(£,C)nu(£,C) (with ^ even) 

T)e (^> 4) : so(2£,M) . 

As £ > 3, the isomorphisms (see, e.g., Thm. X.3.12 in 
[g^l) tti ~ bi ~ Ci and b2 ~ C2 or the isolated semi-simple 
case c)2 ~ ai © Qi are of no concern. The odd-dimensional 



type bi is ruled out as i is an irreducible simple subal- 
gebra to su(2") thus leaving the types q and dg as the 
only classical alternatives to ag. Note that option (3) is 
impossible by the promise of spin-^ representations: e.g., 
in su(4) there is an irreducible spin-| representation of 
su(2) as a proper subalgebra su(2)^-^3^2 £ ^"(4), but a 
change in spin quantum number j and number of qubits 
n is unphysical and ruled out by premiss. So if J is a 
classical proper simple subalgebra to su(2"), then it is 
conjugated to so(2") or to sp(2"~^), in which case the 
{iH^} would be either jointly conjugate (1) to a set of 
skew-symmetric or (2) to a set of symplectic Hamiltoni- 
ans as stated in the Lemma. 

The exceptional Lie algebras tg,^e,Qi can — by 
dimensionality — be excluded to count for irreducible n 
spin-i representations of simple subalgebras to some 
su(2") in multi-qubit systems, as, e.g., in the simplest 
representations [Ml (see also [6^) one has 

02 C S07; f4 c 5026; C sby; er c spsg; es c SO248 ■ 

This ends the (somewhat preliminary) proof of the 
Lemma and the Corollary. ■ 

Remark 1 Finally suffice it to add to the above sketch 
that — with the single exception of Q2 in an example of di- 
mension N =7 which again is of no concern in the qubit 
systems here — they fail to generate groups acting transi- 
tively on the sphere or on M.^ \ {0} as has been shown 
in fsdi] building upon more recent results in fd^] to fill 
earlier work 163, Lfia/- 



Fortunately, in practice the underlying simple Lie al- 
gebra may already be evident from the context: consider 
for instance a typical spin coupling Hamiltonian 

Hs = y^^Ajk{(Tjx'E)(Tkx)+Bjk{(Jjy®aky) + Cjk{(Jjz®Crkz) 

(22) 

or a common quadratic Hamiltonian 

n 

Hb(F) = X! ^jk^'jak + B.jkaja\ + Cjkalak + Djka]a\ , 

(23) 

where the aj and a] denote the respective annihi- 
lation and creation operators satisfying the canonical 
(anti) commutation relations 

bosons: [aj, a^] — Sj^k and [aj,afe] =0 (24) 
fermions: {oj, a^} ~ Sj,k and {aj,ak} =0.(25) 

Therefore we get the natural representations 

spin systems : : {ii/^} G su(£ -I- 1) (26) 

bosonic systems : ce : {iH,^} G 5p{£) (27) 
fermionic systems: Of : {iH,y} £ so(2£) . (28) 
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(®-®-©-(D) • • •)^. 

Figure 7: Example of a periodic spin chain of Heisenberg or 
Ising coupling type in which the spin-type pattern of type- 
wise jointly controllable spins is repeated n' times. 

Thus the ambiguity that arises when given t with the 
only promise it is an irreducible unitary representation of 
some simple compact Lie algebra i C su(2") has its origin 
in physics: the unitary representation as a (pseudo)spin 
system could simulate or encode by conjugation 

1. a true spin system, 

2. a bosonic system, or 

3. a fermionic system. 

This ambiguity is removed as soon as the (pseudo)spin 
Hamiltonians {iH^} can be checked for either allowing 
or for safely excluding a bosonic or a fermionic repre- 
sentation. If not, a purely numerical assessment is more 
tedious. 

Check for Symplectic Representation fsdl] : Determine 
all solutions S := {J;} to the coupled system J* + J = 
with JH^ + HlJ = for all {H^}. If 5 = {0} there 
is no equivalent symplectic representation of {iH^} C I. 
However, if there exists a non-zero element in S that 
is non-singular, then t is a subalgebra to su(A^) that is 
conjugate to sp(-/V/2). 

Check for Skew- Symmetric Representation: Make use 
of the wee-notation [H^l, where the product of matrices 
AXB becomes (B* ® A) vec {X) with vec {X) being the 
column vector of all columns in the matrix X. Check 
whether there is a joint unitary matrix V G U{N) such 
that {V^V (g) yty) vec {Hi) = - vec {H^) for all v. 

Example 3: Lack of Symmetry 

Consider spin chains of Ising or Heisenberg coupling 
type with repeated patterns of spin types than can only 
be controlled typewise. Let the patterns be arranged such 
that no 5*2 mirror symmetry occurs. An example is illus- 
trated in Fig. [71 Then all instances of non-symmetric 
spin chains (and other connected coupling topologies) 
numerically addressed thus far have clearly reproduced 
{ = su(2") by Algorithm 1 once an independent alge- 
braic analysis had ensured that the commutant to the 
generating set of drift and control Hamiltonians {H^}' is 
trivial. 



II. OBSERVABILITY 

Controllable systems are closely related to 'observable' 
ones, as has been known for a long time in classical sys- 
tems and control theory, see, e.g., [tOI- In the context 



of quantum dynamics, the interrelation has been sum- 
marised in the textbook by D'Alessandro j5§|. There the 
notion of observability is in a direct sense. Consider an 
observed bilinear control system (henceforth termed S) 

p{t) = -i[Hu, p{t)] where p(0) = po (29) 

with the measured output by the detection observable C 

y(t) tr{Ctp(i)} • (30) 

Recall that we specified the reachable set in terms of 
subgroup orbits generated by the dynamic Lie algebra, 
so 

pit) e O-Kipo) ^ {K poK-'\ K e-K = Clip t} . (31) 

Then two states represented as density operators pAo and 
PBO are said to be indistinguishable by the observable C, 
if at any time t and under any control u{t) they give the 
same expectation values {C){t) over the entire reachable 
sets, i.e. if 

tT{C^OK{pA0)}^tv{C^OK{pB0)} ■ (32) 

Thus the notion of indistinguishability is closely related 
to the topic of the C -numerical range [tiI . It^ I being de- 
fined as 

W{C,A) ■.= tr{{C''UAU^)\U e SU{N)} . (33) 

In particular, the set of expectation values {C){t) ex- 
actly coincides with what we introduced as the relative 
C -numerical range [t^. 

Wk{C,A) := tr{{C''KAK'<) \K eKC SU{N)} (34) 

with respect to the compact connected subgroup K gen- 
erated by the dynamic Lie algebra ? of the controlled 
system. Therefore one obtains the following: 

Corollary II. 1 Two states pAo and pBo o,re indistin- 
guishable by an observable C if their relative C -numerical 
ranges with respect to the (compact connected) dynamic 
Lie group K := expt C SU{N) coincide 

W^{C,pao) = W-k{C,pbg) ■ (35) 

Now a controlled quantum system is called observable 
whenever two states are indistinguishable if and only if 
they are equal. In Ref. [H^ D'Alessandro showed that 
fully controllable quantum systems are also observable 
{ibid., Proposition 4.2.8). This is the quantum adapta- 
tion of its famous and even stronger forefather in classical 
systems and control theory, where controllability and ob- 
servability have long been established to be in a duality 
relation: i.e. a classical system is controllable if and only 
if it is observable. For details see, e.g., (70l |. 
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A. Constraints by Symmetry 

Here we resume some of the arguments of [Hsl . [75| and 
relate them to symmetries of quantum systems: to this 
end turn the observables C G f)er(iV) into their traceless 
forms C e t)ero(iV) 

C:=C-'^-1n . (36) 

Then each observable comes with an observability space 
with respect to the dynamic Lie algebra t, which in anal- 
ogy to [73| can be written as 

oo 

Vc ■■= ad^j spanjiC} , (37) 

3=0 

where adj operating on C is defined by the linear span 
of j times repeated Lie brackets of elements h e i with 
iC. Due to the Cayley-Hamilton theorem, the infinite 
sum terminates at some finite value. Vc is a subspace of 
su{N), but not necessarily a subalgebra. The observabil- 
ity space (and its orthocomplement Vq of Vc in u{N)) 
both form ideals under the dynamic Lie algebra 

%Vc]CVc and [t,V^]CV^ . (38) 

Moreover in ref. [5§| it was also shown that two states 
PAo and pbo are indistiguishable by C if and only if their 
difference lies in the orthocomplement {pAo — PBo) S i^c- 
This appears very natural, as the difference between dis- 
tinguishable states has to be observable. So the con- 
trolled system S of Eqn. l(29|) with output as in Eqn. l(30|) 
is observable by C if 

Vc '= su(iV) . (39) 

By similar arguments, one may extend the controlled 
bihnear system to be observed by several non-trivial de- 
tection operators C := {Ci, C2, ■ ■ ■ , C*, }, see [t^ for de- 
tails. Since the respective observability space Vc is again 
an ideal under the dynamic Lie algebra t, and full control- 
lability entails ? = su(A^), then t is a simple Lie algebra 
by construction and has only trivial ideals, su(A^) itself 
and the identity — which correspond to Vc and Vq, re- 
spectively. Therefore a fully controllable system is also 
observable [59| . 

Recall that the commutant i' to the dynamic Lie al- 
gebra t constitutes an ideal to i. As in Sec. I, non- 
trivial symmetries preclude that the dynamic Lie alge- 
bra is simple. An immediate consequence is that non- 
trivial symmetries s G t' not only exclude full control- 
lability but also restrict observability, because the or- 
thocomplement Vq- becomes non trivial. This can be 
made more precise by the following natural definition: 
Let S be again a bilinear quantum control system with 
output observed by a given set of detection operators 
C := {Ci, C2, . . . , Cr} that shall be mutually orthogo- 
nal with respect to the Hilbert-Schmidt scalar product. 



Let V(E, t,C) be the associated observability space un- 
der the dynamics of E governed by the Lie algebra t. 
Then the system is called observable by C on the sub- 
space TZ C f)er(iV) if TZ := spaujj {pi, p2, ■ ■ ■ , Pi}, where 
{pi, P2, • ■ • , Pi} is a maximal set of linearly independent 
states such that for all pairs {p„,p^) G TZ x TZ 

(p.-PM)^iV^(E,e,C) . (40) 

Clearly, the observable subspace TZ grows with the num- 
ber of (experimentally) available detection operators C, 
or more precisely, with the number of disjoint orbits un- 
der the dynamic group K, i.e. C'k(C'i/) with G C. Now 
if the observables C share the same symmetry group as 
the dynamic system and hence its algebra t, then the 
observability space must be a subspace to t. 

Therefore in actual experiments we have the following 

Proposition II. 1 Let Y, be a bilinear quantum control 
system with associated dynamic Lie algebra i. Then an 
operator A S l)cx{N) is observable by the set of detec- 
tion operators C if to the associated observability space 
V(E, t,C) there is a set of preparable linearly indepen- 
dent initial states of the system {pi, p2, ■ ■ ■ , pe} spanning 
TZ so that (i) 3A{t) E Ok{A) : A{t) e TZ and (ii) E is 
observable by C on TZ. 

The consequences for quantum system identification 
are immediate and, fortunately, not limited by the in- 
trinsic symmetry of the system: 

Corollary II. 2 A closed controlled quantum system E 
with dynamic Lie algebra t can be identified by a set of 
orthogonal detection operators C := {Ci, C2, . . . , Cr} in 
the following two scenarios: 

1. if the system Hamiltonian Hq is entirely unknown, 
then it is sufficient that the observability space as- 
sociated to the observables C is V(E,6, C) = su(A^); 

2. if only the structure of the system Hamiltonian Hd 
is known in the sense it is a linear combination 
of orthogonal basis elements, where the knowledge 
about the expansion coefficients is confined to 'zero ' 
or 'non-zero ', then it suffices all the basis elements 
with non-zero coefficients are observable by C. 

In a sense, scenario (1) mirrors full controllability, 
whereas scenario (2) refiects 'task controllability' with 
respect to all non-zero Hamiltonian components. — In 
relation to the second scenario note that the dynamic 
Lie algebra t is itself independent of the actual coupHng 
strengths as long as they are sizeably different from zero 
(and do not decide whether a symmetry is introduced or 
is broken). 

Up to now, we have treated observability in one step. 
However, similar arguments generalise naturally to sev- 
eral measurements at k different times [t^I and thus form 
a basis of quantum process tomography. For instance, 
then the notion of being indistinguishable is replaced by 
being indistinguishable in k time steps. 
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Moreover it is to be anticipated that repeated mea- 
surements will play a further role in indirect Hamiltonian 
identification as introduced recently [tH] to determine the 
parameters of a system Hamiltonian by tracing time evo- 
lutions measured on just a subset of gateway qubits. 

B. Constraints by Relaxation 

Consider a controlled bihnear system Ej\f with Marko- 
vian relaxation expressed in GKS-Lindblad form 

p = -iada„(p)-ri(p) , (41) 

where the output shall be measured by the observable C 

y{t)^tr{C^p{t)} . (42) 

Controllability and observability issues in those type of 
controlled open quantum system are more intricate as 
is evident even in the simplest case: Let the iad^r^ be 
skew-Hermitian, while Tl shall simplify to be Hermitian 
(which generically need not hold). The usual Cartan de- 
composition of 0[(A^^,C) := 4 ® p into skew-Hermitian 
matrices (t) and Hermitian matrices (p) with its commu- 
tator relations C t, [p,p] C 6, [6, p] C p illustrates 
that double commutators of the form 

[[iadff^.,ri],[iadjf,,ri]] (43) 

may generate new t-directions in the dynamic Lie algebra 
that by relaxation turns into a subalgebras C a[(A^^ C ). 

As we have recently illustrated on a general scale |77l |. 
controlled open systems thus fail to comply with the stan- 
dard notions of controllability, see also [60, [ill- If, in the 
absence of relaxation, the Hamiltonian system is fully 
controllable, one finds 

(ii/d, iHj I J - 1, ... , m)ue = 5u{N) , (44) 

while master equations like Eqn. (|4T|) generically give a 
dynamic Lie algebra 

Sopcn (iadHd +rL,iadHj | j = 1,2, . . . ,m)Lie 

(45) 

- 0[(fiero(7V)) , 

cf. [sE, [73, [li, [83|. Then the dynamics of the entire 
open system takes the form of a contraction semigroup 
contained in GL{^zXq{N))\ the relaxative part interferes 
with the coherent Hamiltonian part generating new di- 
rections in the Lie algebra, where the geometry of the 
interplay determines the set of explored states. 

Hence we introduced weaker controllability concepts in 
open systems (tJ: fix the subalgebras generated by the 
control terms as 

«c := (ii/i,...,ii/™)Lic (46) 
and the one extended by the Hamiltionian drift as 

id:= {iHd,iHi,...,iH^)uc ■ (47) 
Then an open Markovian quantum system Ea/ I|4T|) is 



(a) Hamiltonian- controllable (h), if ?c = su{N) and no 
bounds on the control amplitutes Uj, j = 1, . . . ,m 
are imposed (idealised instantanous reachability); 

(b) weakly Hamiltonian- controllable (wh), provided 
td=5u{N) and Tl 1 =7! with 7 > 0. 

In general, an open quantum system that is fully con- 
trollable in the absence of relaxation need not be WH- 
controllable when including relaxation. 
Now with regard to observability we obtain 

Corollary II. 3 Let be a bilinear controlled open 
Markovian quantum system following Eqns. l(4T|) and 
(|42|) . Assume that in the absence of relaxation, the sys- 
tem E is fully controllable and thus comes with the asso- 
ciated observability space V(I], t, C) — su{N). Then 

1. if it is E-controllable, the open system has the same 
observability space V(T,M,i,C) = V(S, J, C); 

2. if it is wu- controllable, the open system has an ob- 
servability space that is the contraction of su{N) in 

3. generically, the open system has an observability 
space that is a subset in g[(f)ero(iV)) growing with 
time after which the measurements are taken. 

We recently gave an illustration of the corresponding 
reachable set in generic open Markovian systems [73| 
showing how intricate it is to specify reachable sets in 
generic open systems. Hence addressing observability 
spaces in open systems is an equally subtle enterprise 
with lots of open research problems. 



III. CONCLUSIONS AND OUTLOOK 

We have treated controllability and observability in 
a unified Lie-algebraic framework incorporating con- 
straints by symmetry for closed systems and constraints 
by relaxation in open ones. In particular, the dynamic 
system Lie algebra allows for specifying the reachability 
sets of closed systems explicitly. In symmetric quantum 
hardware setups the feasible tasks in quantum simula- 
tion or quantum gate synthesis can be made precise thus 
giving valuable guidelines for quantum hardware design 
matched to solve given problems. 

Likewise, the dynamic system Lie algebra determines 
those Hilbert subspaces that can be directly observed 
with a given set of detection operators. We anticipate 
that the approach will prove widely useful in Hamiltonian 
system identification as well as in process tomography. 

Using the Lie cone as the open-system generalisation 
of the dynamic system Lie algebra in closed ones seems 
a logical consequence, yet in practice it is by no means 
straightforward. Thus in [t^I we gave a sketch how to 
approach it from the inside as well as from the outside. 
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